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Abstract —  This  paper  focuses  on  computing  and  interpreting 
the  significance  of  the  observability  Jacobian  singularities. 
Previous  efforts  have  focused  on  developing  an  observability 
index  and  determining  the  effects  of  model/measurement  changes 
to  the  observability  formulation.  Using  this  index,  the 
singularities  of  the  Jacobian  are  extracted,  and  preliminary 
results  indicate  that  there  is  a  duality  between  the  loss  of  stability 
(whether  it  be  phase  or  voltage  stability),  and  the  loss  of 
observability.  Thus,  providing  a  single  dual-use  metric  capable 
of  forecasting  both  stability  and  observability  characteristics  of  a 
system.  The  results  are  illustrated  using  a  sample  3-bus  system. 

Index  Terms — observability,  power  systems,  system  index, 
differential  algebraic  equations. 

I.  Introduction 

The  safe,  reliable  and  economic  operation  of  modern  power 
systems  rely  on  increased  levels  of  system  automation, 
power  electronic  equipment,  and  system  controllers  capable  of 
monitoring  and  controlling  the  system  through  various 
operating  points. 

One  particular  family  of  power  systems  where  the  effects 
of  system  automation,  power  converters  and  rapid  system 
reconfiguration  are  more  profound  are  naval  shipboard  power 
systems.  The  inherent  nature  of  the  shipboard  systems, 
similar  to  island  topologies  in  terrestrial  systems,  has 
magnified  the  need  for  real-time  system  control  and  therefore 
the  development  of  dynamic  state  estimators  and  nonlinear 
observers.  Shipboard  power  systems  are  unique  in  their 
composition,  and  differ  from  terrestrial  systems  in  many  ways 
[1]: 

•  Single  dynamic  loads  consuming  a  large  portion  of  the 
generation  capacity 

•  Time  constants  of  the  generators  are  closer  to  the 
electrical  time  constants  than  in  terrestrial  systems 
•  System  reconfiguration  routines  are  much  faster  and  a 
large  portion  of  the  system  load  is  sensitive  to  power 
interruptions.  A  small  interruption  of  approximately  100 
msec,  can  cause  subsystems  to  shut-down. 

•  Redundancy  is  incorporated  in  all  aspects  of  system  design 


to  ensure  survivability 

The  power  demand  on  a  ship  has  evolved  from  kilowatts 
the  tens  of  megawatts,  and  is  constantly  increasing  with  new 
ship  designs  and  concepts.  One  particular  concept  affecting 
directly  the  shipboard  power  system  is  the  vision  for  an  All 
Electric  Ship  [2],  which  signifies  the  move  of  the  shipboard 
power  system  from  an  electromagnetic/electromechanical 
basis  to  a  power  electronic  basis.  Some  of  the  design  aspects 
involving  the  All  Electric  Ship  concept  are: 

•  Electric  Drive 

•  Integrated  Power  System 

•  Zonal  Distribution  System 

•  Reduced  Manning 

The  enabling  technology  in  all  of  the  above  concepts  is 
intelligent  automation  and  control.  With  the  increase  in  the 
automation  of  system  functions  within  the  shipboard  power 
system,  the  dependence  on  both  executive  and  distributed 
controllers  to  maintain  the  performance  of  the  power  system  at 
acceptable  levels  is  becoming  more  apparent.  Traditional 
control  approaches  based  on  linearized  models  of  the  system 
generation/distribution  capacity  will  not  be  sufficient.  The 
need  of  incorporating  the  nonlinear  dynamics  of  the  shipboard 
power  system  in  the  control  methodology  is  the  solution  to 
this  problem,  and  one  proposed  solution  is  the  use  of 
nonlinear  observer  based  controllers  and  dynamic  state 
estimators.  Before  these  solutions  are  designed  however,  the 
concept  of  nonlinear  observability  is  first  investigated. 
Providing  an  implementation  of  the  observability  formulation 
for  power  systems  modeled  as  Differential  Algebraic 
Equations  is  the  starting  point  for  the  development  of 
nonlinear  controllers,  dynamic  state  estimators  and 
performance  index  measures. 

The  authors  have  investigated  the  problem  of  system 
observability  for  nonlinear  power  systems  modeled  as  DAEs 
[3-6],  and  this  work  emphasizes  on  the  possible  duality 
between  the  singularities  of  the  observability  Jacobian  and  the 
known  stability  characteristics  of  the  system.  The  preliminary 
results  indicate  that  a  system  becomes  unobservable  at  the 
maximum  load  point. 
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II.  Problem  formulation 

The  majority  of  the  system  models  used  in  this  work  and 
presented  in  this  section  are  adopted  from  [4]  and  repeated 
here  for  the  sake  of  continuity.  The  Differential  Algebraic 
Model  of  a  power  system  consisting  of  n  generators  and  m 
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static  load  buses  is  given  as: 

x  =  f(x,u,N )  / : R2(',+m“1)  R2("_1) 

0  =  g(x,u,N )  g;]R2("+m-1)  ^R2(m)  (1) 

y  =  h(x,N )  h :  R2("+ra_1)  IT 

where  /(.)  is  the  set  of  differential  equations  related  to  the 
dynamics  of  the  generation  equipment,  g(.)  is  the  set  of 
nonlinear  algebraic  equations  related  to  the  static  load  buses 
of  the  system,  h(.)  are  the  measurement  equations,  u  are  the 
system  inputs,  N  describes  the  characteristics  of  the  individual 
power  lines  in  the  system,  v  is  the  system  vector  containing 
both  the  dynamic  states  and  static  variables  of  the  system,  and 
y  are  the  available  measurements.  More  information  on  DAE 
power  system  models  and  their  related  characteristics  can  be 
found  in  [7]. 

By  reformulating  (1)  in  a  more  condensed  version: 

F(x,x,N)  =  Bu  (2) 

y  =  h(x,N) 


the  observability  Jacobian  {J0)  of  the  system  is  given  by  [8]: 
~G  G  G 

X  X  V 

7/  H{  H , 


(3) 


where 


F(x,  x,N) 

Fm  :  Fx (x, x, N)Jc  +  F. (x,  x, N)x 


F(s)  :(F(x,x,N))M 


h(x,N) 

y 

H  = 

hm  :  hx(x,N)x 

= 

ym 

_h(r)  :h(x,N)(r\ 

_y(r) _ 

(4) 


and  s  is  the  system  differentiation  index,  r  is  the  measurement 
differentiation  index,  and  w  are  the  higher  derivatives  of  x 
(2....w).  Based  on  this  Jacobian,  the  system  is  said  to  be 
smoothly  observable  if  [8]: 


1. 


rank(J0 )  =  n  +  rank 


G, 

Hx 


Gw 

Hw 


2.  rank(J0)  is  constant  rank  on  S 


(5) 


S:{0  =  g(x,u,N)} 


A.  Power  System  Model 

In  the  adopted  structure  preserving  DAE  model  of  the 
power  system,  the  generators  are  modeled  using  the  rotational 
dynamics  provided  by  the  classical  swing  equation,  and  the 
loads  are  modeled  as  constant  PQ  loads  (considered  power 
distribution  centers),  forming  the  system  equations: 

f1  :Si-ct)i=  0 

/"  :  a,  +  A (^ )  +  JL  P  (s,  9,V,N)  =  ^-  (6) 

Mf  '  M{ 

fm  :Pj -Pj(S,0,V,N)  =  O 
fw  :Q]-Qj{S,e,V,N)  =  Q 


N  -Ycd  =  Gcd  +  jBcd 


i  =  2...n  generators;  j=n  +  l...m  loads; 
c,  d  =  +  m 


where  Sh  Mh  Dt  and  PMi  are  the  internal  angle,  inertia, 
damping  and  constant  mechanical  power  of  the  ith  generator 
respectively,  and  Oj  and  V}  are  the  jth  load  bus  angle  and 
voltage  respectively.  The  desired  power  demand  on  the  jth 
load  bus  is  denoted  S*=P*+jQ*.  The  system  has  the  form  of 

(1) ,  with  f(.)  in  (1)  comprised  of  /  and  f  and  g(.)  in  (1) 
comprised  of  fn  and  /y,  and  with  respect  to  the  form  given  in 

(2) : 


F(x,a)=[/7/rnr 

x  =  [Si,col,0J,VJ],u  =  ^- 


(7) 


The  available  measurements  are  all  assumed  as  real  power 
branch  flows,  given  as: 

h(x,N)  =  {Pcd(x,N)}  (8) 


This  system  can  be  further  modified  to  incorporate 
generator  exciter  dynamics  or  dynamics  in  the  system  loads. 
The  observability  formulation  for  these  variations  was 
discussed  in  [3-6],  and  the  addition  of  other  types  of  generator 
and  load  models  (voltage  dependent  loads  for  example)  will 
impact  the  construction  of  the  observability  Jacobian  only. 
The  methodology  for  obtaining  the  Jacobian,  the  observability 
conditions  and  the  degree  of  observability  measure  is 
consistent.  Also,  it  should  be  noted  that  in  (6)  generator  bus  1 
was  arbitrarily  selected  as  the  reference  bus,  yielding  SgSrS] 
and  Of  6}- 8]. 


B.  Observability  Metric 


To  quantify  the  observability  of  the  system  during  a 
variation,  a  measure  of  observability  is  introduced,  defined  as: 

A- (4)  l  =  l...L  (9) 


n  =- 


K**  (Jo) 

where  are  the  singular  values  of  the  observability  Jacobian 
obtained  via  a  singular  value  decomposition  at  each  point  of 
interest  /.  The  maximum  number  of  L,  is  dictated  by  the 
convergence  of  the  numerical  solvers  used  for  either  static 
studies  (points  along  PV  Curve),  or  dynamic  studies  (length  of 
time  window  in  a  time  domain  simulation)  . 

The  degree  of  observability  y  is  defined  as: 

y=  maxjV]  l  =  l...L  (10) 


Based  on  the  formulation  of  this  measure,  a  smaller 
numerical  value  of  /is  preferred  -  a  higher  number  indicates 
that  the  observability  Jacobian  is  becoming  singular. 


III.  Approach 

To  determine  all  the  singularities  of  the  observability 
Jacobian,  a  symbolic  computation  approach  is  used.  Also, 
time  domain  simulations  are  used  to  track  the  observability  of 
the  system  along  the  system  state  trajectories,  as  these 
unobservable  regions  are  approached.  A  brief  description  of 
both  the  symbolic  and  numerical  approaches  is  provided. 
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A.  Symbolic  Approach 

The  symbolic  approach  examines  the  observability  criteria 
as  a  function  of  the  system  state  space.  In  particular,  the 
examination  is  focused  on  finding  where  in  the  state-space  the 
observability  Jacobian  undergoes  rank  deficiencies.  In 
essence,  by  examining  where  observability  is  lost,  all  other 
points  in  the  state-space  are  assumed  to  be  observable.  This  is 
assumed,  rather  than  guaranteed,  because  of  the  performance 
of  the  symbolic  tools  in  evaluating  rank  deficiencies.  Namely 
their  performance  in  locating  higher  order  degeneracies  in  the 
Jacobian  is  of  question.  For  the  particular  system  examined, 
the  symbolic  tools  were  able  to  capture  higher  order 
degeneracies.  This  technique  is  powerful  in  the  sense  that  it 
provides  the  global  solutions  where  the  system  is 
unobservable,  however  it  is  limited  by  the  computational 
capability  of  the  symbolic  solvers,  and  hence,  is  limited  to  the 
size  of  the  system  that  can  be  evaluated.  The  symbolic 
approach  process  is  outlined  in  Figure  1 . 


Fig  1.  Investigating  Observability  using  Symbolic  Tools 


Fig  2.  Flow  Diagram  of  the  Numerical  Approaches  in 
Determining  Observability 


B.  Numerical  Approach 

The  numerical  approach  consists  of  two  separate  cases: 

•  Static  Cases  -  Evaluating  observability  for  operating 
points  along  system  PV  Curves 

•  Dynamic  Cases  -  Evaluating  the  observability 
criterion  along  trajectories  of  the  system  states 

The  general  approach  in  evaluating  these  criteria  and 
determining  the  appropriate  s  and  r  to  render  the  system 
observable  is  outlined  in  Figure  2. 

For  the  static  cases,  the  solution  of  the  system  PV  curves  is 
accomplished  in  an  iterative  fashion  using  the  Newton- 
Raphson  (NR)  and  Newton-Raphson-Seydel  (NRS) 
algorithms  [9].  When  the  NR  algorithm  fails,  the  NRS 
algorithm  is  used  to  obtain  the  data  points  near  and  around  the 
nose  point  of  the  PV  curve,  as  illustrated  in  Figure  3. 

For  the  numerical  cases,  the  ode  15s  DAE  solver  was  used 
in  the  Matlab/Simulink  environment  to  generate  the  time 
domain  simulations  of  interest. 


Fig  3.  Typical  Power  System  PV  Curve 

IV.  Simulation  Results 
A.  Sample  System 

To  illustrate  the  symbolic  approach  for  obtaining  the 
unobservable  regions,  a  3 -bus  power  system  depicted  in 
Figure  4,  with  corresponding  system  parameters  in  Table  1, 
was  used.  There  are  two  sets  of  results  -  considering  two 
different  measurements,  Pi2  and  P32.  In  both  cases,  r=s=l , 
and  the  singularity  of  the  Jacobian  is  investigated  by  solving 
the  following  relationship  symbolically: 

det  ( )  =  0 

For  the  selected  system,  the  observability  Jacobian  has  the 
general  form  provided  in  (1 1),  and  the  states  of  the  system  are 
given  as: 

X  ~  \@2  » »  ^3  >  ^3  ] 
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Fig.  4.  Sample  3-bus  Power  System 
TABLE  I 

System  Parameters  for  Sample  3-bus  Power  System 


Line 

R 

X 

D=  1 

1-2 

.0194 

.0592 

M=  1 

2-3 

.0470 

.1980 

1-3 

.054 

.2230 

Bus 

pinitial 

Qinitial 

Vi=1.06 

2 

.1830 

- 

V2=1.04 

2 

-.9420 

.0440 

B.  Symbolic  Results  -  Loss  of  Observability 

By  computing  the  determinant  symbolically  and  equating  it 
to  zero,  there  are  6  solutions  for  x  which  yield  rank 
deficiencies  (Table  II).  In  analyzing  the  solutions,  and 
keeping  in  mind  the  significance  of  the  respective  system  state 
assignments,  the  first  five  solutions  are  not  considered 
physically  realizable.  For  Solution  1,  the  bus  voltage  at  the 
load  bus  would  have  to  go  to  zero,  meaning  the  generators  are 
either  offline  or  the  topology  of  the  system  has  changed  (due 
to  a  fault  on  Bus  3,  or  faults  on  both  lines  connecting  Bus  3  to 
the  two  generator  buses).  Similarly,  for  Solutions  2-5,  the 
angle  difference  between  generator  1  and  the  load  bus  is  90 
degrees,  which  would  indicate  a  phase  instability  in  the 
system.  Therefore,  for  the  first  five  solutions,  one  can  assume 
that  the  system  will  be  smoothly  observable  for  all 
perturbations  that  do  not  cause  these  specific  system 
instabilities,  and  that  one  measurement  (P32)  is  sufficient  for 
the  system  to  be  smoothly  observable.  The  only  solution  of 
interest  now  is  Solution  6,  which  represents  the  load  bus 
voltage  as  a  function  of  the  two  bus  phases  in  the  system. 
Similarly  for  the  case  of  a  measurement  P12  (Table  III),  the 


first  2  solutions  represent  either  a  fault  on  bus  3,  or  a  phase 
instability,  and  the  third  solution  provides  the  bus  voltage  as  a 
function  of  the  bus  phases. 

TABLE  II 

System  Observability  is  Lost  and  the  Corresponding 
Rank  of  the  Observability  Jacobian  -  P32 


Solution 

Xl 

x2 

x3 

X4 

Rank 

l 

Xi 

x2 

X3 

0 

6 

2 

0 

X2 

71/2 

X4 

8 

3 

0 

X2 

-71/2 

X4 

8 

4 

71 

X2 

71/2 

X4 

8 

5 

n 

X2 

-71/2 

X4 

8 

6 

Xi 

X2 

X3 

f(Xi,  X3) 

8 

TABLE  III 

System  Observability  is  Lost  and  the  Corresponding 
Rank  of  the  Observability  Jacobian  -  P12 


Solution 

Xl 

x2 

x3 

x4 

Rank 

l 

Xi 

X2 

X3 

0 

6 

2 

71/2 

X2 

X3 

x4 

8 

3 

Xl 

X2 

X3 

f(Xl,  X3) 

8 

Examining  Solution  6  of  Table  II  and  Solution  3  of  Table 
III,  and  similar  solutions  for  a  lossless  system,  an 
unobservable  region  can  be  constructed,  given  a  range  for  the 
bus  angles  is  specified.  To  examine  the  physical 
interpretation  of  these  solutions,  the  initial  operating  point  in 
Table  I  is  selected  with: 

jc  =  [-.0869,0, -.2220, 1.01] 

and  an  arbitrary  PV  curve  is  constructed  (Figure  5).  Using  the 
phase  quantities  dictated  by  the  PV  curve,  the  relationship  of 
the  load  bus  voltage  to  the  phase  difference  between  bus  2  and 
3  is  examined  and  compared  to  the  results  for  the  load  bus 
voltage  obtained  using  the  PV  curve.  The  phase  difference 
was  used  as  a  parameter  to  condense  the  problem  into  2-D. 

Examining  Figure  6,  the  solution  for  the  load  bus  voltage 
magnitude  intersects  the  unobservable  solutions  obtained  from 
Tables  II  and  III,  and  using  the  information  obtained  by  the 
PV  curve,  this  intersection  corresponds  to  the  maximum  load 
point  of  the  system  (magnified  in  Figure  7).  This  result  is 
consistent  with  all  the  obtained  solutions,  and  is  invariant  to 
the  selected  measurement,  and  to  the  implication  of  a  lossless 
or  lossy  system.  Due  to  the  symmetry  of  the  system  topology 
used,  the  results  for  P13  are  the  same  as  those  for  P32  and  are 
not  included. 


PV  Curve  for  Sample  3-Bus  System 


Fig  5.  Arbitrary  PV  Curve  of  3-Bus  System 
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Relationship  Between  System  PV  Curve  and  Unobservable  Regions 


Fig.  6  Relationship  Between  System  PV  Curve  and 
Unobservable  Solutions 


Relationship  Between  System  PV  Curve  and  Unobservable  Regions 


Theta2-Theta3  (rad) 

Fig.  7  Relationship  Between  System  PV  Curve  and 
Unobservable  Solutions  -  Intersection  Magnified 


along  the  trajectories  of  the  system  states.  It  is  expected  that 
the  observability  index  for  the  trajectories  associated  with 
point  C  will  be  higher.  By  repeating  this  process  for  the  entire 
set  of  points  along  the  PV  curve,  a  threshold  on  the 
observability  index  was  obtained  (j*=  600),  and  the  individual 
/is  provided  in  Table  V,  for  the  perturbation  level  of  40%  on 
the  generator  phase. 


TABLE  IV 

Arbitrary  Points  along  PV  curve  chosen  for  time 

DOMAIN  SIMULATIONS 


Point 

Xi 

x3 

x4 

Pi 

Pa 

q3 

A 

-.0076 

-.1167 

1.0218 

.3830 

-1.1420 

.0440 

B 

.0461 

-.5918 

.7196 

3.4006 

-4.1596 

.0440 

C 

.0047 

-.6274 

.6911 

3.4302 

-4.1892 

.0440 

Point  A 


Time  (sec.) 


Therefore,  for  the  selected  system,  a  duality  exists  between 
the  loss  of  observability  and  the  voltage  stability  problem.  At 
the  maximum  load  point,  the  system  is  unstable  and 
unobservable.  This  is  also  present  with  cases  of  phase 
instability,  represented  by  the  solutions  in  Table  II  and  III, 
where  large  phase  differences  exist  between  the  generator  and 
load  bus. 

C.  Dynamic  Simulations 

Because  of  the  limitations  with  symbolic  tools,  analyzing 
larger  systems  are  not  possible  without  the  introduction  of 
numeric  methods.  However,  unlike  the  symbolic  approach 
which  yielded  the  entire  family  of  unobservable  solutions,  the 
numerical  approaches  are  constrained  to  only  approach  a  point 
of  interest.  To  illustrate  that  the  numerical  results  are 
consistent  with  the  answers  from  the  symbolic  computations, 
three  arbitrary  points  are  selected  along  the  PV  curve  (points 
A,  B  and  C  in  Figure  5,  given  in  Table  IV).  A  40% 
perturbation  is  applied  to  the  generator  phase  at  each  point, 
and  the  resulting  time  domain  waveforms  are  extracted 
(waveforms  provided  in  Appendix).  These  time  domain 
waveforms  are  then  used  to  analyze  the  observability  index 


Fig.  8.  Observability  Measure  along  State  Trajectories  for 
Points  A,  B  and  C. 

TABLE  V 

Observability  Index  for  Each  Set  of  Trajectories 
Initiated  from  Points  A,  B  and  C 


Point 

Y 

A 

101.6915 

B 

124.4947 

C 

265.4845 

As  expected,  the  trajectories  produced  from  the  points 
closer  to  the  nose  point  produced  larger  index  values, 
indicating  that  the  system  is  becoming  less  observable,  and 
they  are  consistent  with  the  results  obtained  using  the 
symbolic  approach.  However,  at  no  point  does  the  system 
become  unobservable  in  a  numerical  sense.  The  observability 
criteria  of  Equation  5  are  satisfied  for  all  the  points  along  the 
PV  curve,  and  system  trajectories  for  this  system.  A  closer 
examination  indicates  that  a  quantitative  change  is  occurring 
as  the  system  is  approaching  the  maximum  load  point.  This  is 
evident  when  examining  the  minimum  singular  value  of  the 
observability  Jacobian  for  each  point  along  the  PV  curve.  As 


Trajectories  After  Perturbation  from  Point  B 


illustrated  in  Figure  9,  the  value  approaches  zero  at  the 
maximum  load  point.  From  a  numerical  standpoint,  the 
system  is  still  observable,  and  this  problem  is  hence 
transformed  into  a  numerical  sensitivity  issue.  The 
appropriate  selection  of  thresholds  on  y  provides  a  better 
depiction;  however,  this  process  involves  extensive  system 
simulations  to  set  these  thresholds. 


Minimum  Singular  Value  of  Jacobian  Along  PV  Cun© 


Fig  9.  Minimum  Singular  Value  of  Jacobian  Along  PV  Curve 
V.  Conclusions 

Singularities  of  the  observability  Jacobian  correspond  with 
system  singularities  that  are  typically  associated  with  phase 
and  voltage  stability  problems  in  power  systems.  This  duality 
can  be  exploited  by  using  the  observability  index  as  a  system 
metric  capable  of  forecasting  both  the  stability  and 
observability  classifications  of  a  power  system.  The  dual  use 
of  a  single  metric  is  very  encouraging.  In  the  case  of  system 
stability,  since  the  common  points  are  fewer  than  the  entire 
family  of  unstable  points  for  a  system,  this  metric  would  serve 
as  a  necessary  but  not  sufficient  condition. 


Fig  11.  Trajectories  After  Perturbation  from  Point  B 


Fig  12.  Trajectories  After  Perturbation  from  Point  C 
VII.  References 


VI.  Appendix 

Time  domain  simulations  for  the  sample  system  considered 
in  previous  section  are  provided  in  Figures  10-12. 


Trajectories  After  Perturbation  from  Point  A 
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Fig  10.  Trajectories  After  Perturbation  from  Point  A 
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